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ON REPEATED ADJUSTMENTS, AND ON SIGNS OF 
RESIDUALS. 



BY E. L. DE FOREST, M. A., WATER-TOWN, CONN. 

Let any consecutive terms in a series be denoted by 
.... M_j, W—i, U(j, u lt m 2 , . . . &c. 
Suppose them to be adjusted by the formula 

and let the adjusted terms be again adjusted by the formula 

MO = I M 0+ i lK+"-l)+ i 2K+ {( -2)- ( 2 ) 

By the first operation, five consecutive adjusted terms will be 

«L 2 = Zo W -2 + M W -l+ W -3) + Z 2K+ W -4) + *3( M l+ W -5)+*4( M 2+ W --6) 

U ' =l U +?!(«! +U_ 1 ) + l 2 (u 2 +U_- i ) + l 3 (u 3 -\-U_ 3 ) + l 4 (u i +U_ 4 ) 
U \ = Z «1 +Z 1 (M 2 +M ) + ^ 2 ( M 3+ W -l) + ? 3( W 4+ U -2) + ? 4( M 5+ M -3) 
«2 =^0 W 2 + *l( U 3 + W 1 ) + h( U i+ U ) + h( U 5+ U -l) + h( U 6+ U -2)> 

and by the second operation, a readjusted term will be 
Substituting in this the values of u' , u\, &c, as above, we get 



K+«_i) 

(m 2 +m_ 2 ) 



(3) 



u" = [i ? +i 1 (z 1 +? 1 )+i 2 (? 2 +? 2 y 
+ [i ? 1 +i 1 (? +z 2 )+i 2 (? 1 +z 3 ); 

+ [& *2+-&l(*l+*8) + -&aUo+01. . 

+ [i ^3+ i l( Z 2+ ? 4) + - i 2^]( W 3+ W -3) 

+ (-^0^4 + ^1^ +i 2 ^)( W 4+ M -4) 

+ ( i 1^4+ i 2 ? 3)( W 5 + M - 5 ) + - E 2 ? 4K+ M -6> 

It thus appears that the re-adjusted term «' ' is a linear function of the ori- 
ginal terms u , u x , &c, and that the expression (3) may itself be regarded 
as an adjustment formula which replaces and is equivalent to the two form- 
ulas (1) and (2) used successively. Furthermore, if the order of procedure 
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is reversed, and the original series is adjusted, first by (2), and this again by 
(1), it will be found that the re-adjusted term is given by precisely the same 
formula (3). The operations of successive adjustment, then, are of such a 
kind that the order of operation is immaterial. This has already been re- 
marked by me in the Smithsonian Report of 1873, p. 337, though coupled 
with an erroneous inference about the probable error of such an adjustment, 
for correction of which see Analyst for May 1877, p. 84. The fact that 
(3) replaces (1) and (2) without regard to the order of their application is 
one which holds good for any values we please to assign to the coefficients 
l 0) l lt &c, and L , L 1} &c, and also holds good without regard to the na- 
ture of the series u , u 1} &c. The coefficients in (3) are functions of those 
in (1) and (2), and may be computed when these are given. If (1) and (2) 
are adjustment formulas each of which gives an exact result when applied 
to a series of equidistant terms of an order not higher than the third, that 
is, a series whose third differences are constant, then the terms of such a 
series adjusted by (1) will be the same as the original terms, and these re- 
adjusted by (2) will still remain the same, so that (3), which is equivalent 
to (1) and (2) together, will be also an adjustment formula which gives exact 
results when applied to such series. Its coefficients, therefore, will be such 
as to satisfy the two relations (2) given at p. 80 of the Analyst referred to, 
these relations being satisfied, in the case supposed, by the coefficients of the 
component formulas (1) and (2). 

Results similar to the above will be obtained when formulas other than 
(1) and (2), that is, formulas including some other numbers of terms, are 
successively used. We may generalize them as follows. Let the two com- 
ponent formulas be 

m'o =h u o +h ( U l+ U -l)+h («2+«- 2 )+ + l ™ («m + «-m) (4) 

K — ^0 W 0+i 1 («l+M-l) + i 2 ( W 2+ W ~ 2 )+ +■£» K +«-»), (5) 

and for convenience let the second one be the shorter, or n < or = m. A 
comparison of (3) shows that the resultant formula will be 



:={_L l +L 1 {l 1 +l 1 )+L 2 (l 2 +l 2 )+....+L n (l n +l n )■ 

+ [Vl+£l(*0+*2)+-£2(*l+y+--"+-£.&-l+Wl)' 

+ lL l 2 +L 1 (l 1 +l 3 )+L 2 (l a + l 4 )+....+L n {l n _ 2 +l n+2 ) 
+ tL l s +L 1 (l a +l i )+L a {l 1 +l B )+....+LJ! n _ a +l+ s y 



Oo 

(Wl+M-l) 

(u 2 +u_ 2 ) 
(w 3 +«_<,) 



... (6) 
The mode of construction of this expression is easily seen. The coefficients I 
are arranged in columns from the top downward, in the order in which they 
stand, m the series, l m) <„,_!, f m _ 2 , ■ • • 1$, '2? ^i> ^o> 'u ^2> ^3> ■ • • C-u C- 



+ 


L 


U> 




"\ 


+ 


^m- 


-iL) K 


+ w - 


-l) 


+ 


*-'m- 


-2 U K 


+ M_ 


-■?,) 


+ 


L. 


-3<m) ( W S 


+ M_ 


-.0 


+ 


L- 


■ 4 U K 


+ «_ 


-4) 


+ 


<'m- 


- 5 U («£ 


+«- 


- 5 ) 



M7) 
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It must be remembered that all values of I beyond l m are zero, and also 
that wherever the expression gives to I a negative index figure, it must be 
changed to positive, since l_ t =■ l t . 

In this way, by taking m = 4 and n = 2, we find that (6) reduces to (3). 

If a series is adjusted twice by the same formula, then (5) is identical 
with (4), and (6) is somewhat simplified, because we now have 

so that (6) may be written as follows. 

«o = V*h+2(hh+hh+hh+hh+ ■ • 

+ UA+ 2( Joh + hh+hh+hh+ ■ 
+2{l 1 l 2 +l l 3 + l l l 4 +l 2 l 5 + . 

+ UA+Wih + l oh + kh+hh+ ■ 

~T*<'m—l''m\ u 2»i—l T*- (2m— 1)) T 'm'm(. W 2>i~r' M '— 2"J - 

In the successive products ZZ in the horizontal lines of this formula, the se- 
quence of the indices of the first Z, and also of those of the second I, is easily 
seen. Let us apply this to a simple case. The 5-term formula of Table A. 
(Analyst, May 1877, p. 81), in exact numbers is 

"o = M 17u o +12K +u_ 1 )—S(u 2 +u_- 2 )2 (8) 

so that the coefficients I, in decimals of five places are 

l = .48572 l ± = .34286 Z 2 = —.08572. 

Forming the products 

Vo> Vi> ^o^2> n'x> kk) kk) 
and putting m — 2, we get by formula (7) 

u'l = .48572w +.27430(w 1 -f«_ 1 )+.03427(w 2 +w_ 2 ) 

— .05878(w 3 +r{_ 3 )+.00735(w 4 -|-M_ 4 ). (9) 

This is an adjustment formula such that if a given series is adjusted by 
it once, the result is the same as though it had been adjusted by (8) twice. 
Again, if we suppose a series to be adjusted by both (8) and (9) successively, 
the equivalent single formula is found by designating their coefficients by 
L and I respectively, forming all the products of every L by every I, and 
putting m = 4 and % = 2. Formula (6) then gives 
m' "=.41814m +.29305(« 1 +w_ 1 )+.04827(w 2 +u_ 2 )~ .03780(m 3 +m_ 3 ) 
— .01952(« 4 +w_ 4 )+.00756(« 5 +u_ 5 )—.0Q06S(u 6 +m_ 6 ), 

.... (10) 
and this result is denoted by w' " because one adjustment of a series by this 
formula is equivalent to three successive adjustments by the original formula 
(8). In the same way we may combine (10) with (8), (9) or (10), and thus 
get the equivalents of four, five and six applications of (8). The process 



has indeed been extended much farther. It is not so difficult or tedious as 
might be supposed, for it has been found that the coefficients of the extreme 
terms of the formula soon become so small as to disappear, that is, they do 
not affect the first five places of decimals. Thus in (10) the coefficient of 
m 6 is only seen in the fourth place. At the successive steps of the process 
the correctness of the work is readily tested by the condition that for every 
formula obtained we must have 

l +2(l 1 +l 2 +l 3 +&c.) = l, 
within such narrow limits of error as are to be expected from the use of 
logarithms in forming the products. And at each step, before proceeding far- 
ther, it is best to make the formula satisfy this condition exactly, by altering 
one, two or three of the coefficients, if necessary, to the extent of a single 
unit in the fifth decimal place. The various products have been carried to 
five, and sometimes to six, places of decimals. The resulting coefficients 
are shown in the accompanying table I., to three places only. The num- 
bers of successive applications of the original formula (8) are found in the 
upper line of the table, and the corresponding coefficients are in the columns 

Table I. 





1 


2 


3 


4 


5 


6 


8 


12 


16 


32 


h 


.486 


.486 


.418 


.396 


.372 


.354 


.330 


.296 


.274 


.230 


h 


.343 


.274 


.293 


.280 


.277 


.271 


.261 


.246 


.234 


.205 


h 


—.086 


.034 


.048 


.077 


.090 


.102 


.116 


.131 


.139 


.145 


h 




—.059 


—.038 


—.034 


-.023 


-.014 


.002 


.026 


.042 


.072 


h 




.008 


—.019 


—.024 


-.030 


-.032 


-.031 


-.024 


-.015 


.014 


h 






.008 


.000 


-.004 


-.009 


-.017 


-.024 


-.026 


-.017 


h 






—.001 


.004 


.004 


.003 


.000 


-.008 


-.014 


-.022 


h 








—.001 


.001 


•002 


.003 


.002 


-.002 


-.014 


i 










-.001 


-.000 


.001 


.003 


.003 


-.004 


h 


















.002 


.002 


10 




















.003 


hx 




















.001 


r 


.697 


.629 


.596 


.574 


.558 


.544 


.524 


.497 


.479 


.438 


P 


.297 


.102 




.0513 




.0331 




.0160 




.0058 



below. The coefficients in the last column, for instance, are those of an ad- 
justment formula such that if a given series is adjusted by it once, the result 
is the same as if it had been adjusted by the original formula (8) thirty-two 
times in succession, or by (9) sixteen times. This 32-fold formula, if given 
exactly, would include no less than 129 terms, with 65 different coefficients, 
but only 12 coefficients are large enough to appear in the third decimal 
place, and only four more would affect the fourth place. At the bottom of 
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the table are given the ratios of error and of irregularity 



r = 



P W 



The ratio r diminishes slowly as the number of repetitions is increased, while 
p diminishes much faster. For any large number of repetitions, the form 
of the series of values of I approximates somewhat to that of the formulas 
of Table B. (Analyst, p. 82.) For instance, the large coefficients of the 
12-fold and 32-fold formulas of Table I. approach very near to those of the 
11 -term and 15-term formulas of Table B., and their ratios r and p are also 
nearly alike. 

A similar process of repetitions has been applied to the 5-term formula 
of Table B., which in its exact form is 

< = 1 Mlllu + 56(» 1 +«u 1 )—14(» 3 +»_,)]. (11) 

(Smithsonian Report of 1871, p. 333.) The resultant coefficients are given 
in the first half of Table II. The second half contains similar results of 
repetitions of the formula, found at p. 289 of the Report, 



u; 



— ToE 4 ( w o + M i + u -i)— K +«-a)]- 



(12) 



Table II. 





1 


2 


4 


6 


12 


1 


2 


4 


6 


12 


h 
h 

h 

h 

h 

h 

h 
I, 

h 
h 


.570 

.287 

-.072 


.500 
.285 
.001 
-.041 
.005 


.416 
.285 
.059 
-.038 
-.019 
.003 
.002 


.374 

.276 

.089 

—.023 

—.029 

—.004 

.003 

.001 


.310 

.253 

.125 

.015 

—.028 

—.021 

—.004 

.003 

.002 


.400 

.400 

-.100 


.500 
.240 
.080 
-.080 
.010 


.392 

.264 

.101 

-.037 

-.022 

-.008 

.008 

-.002 


.344 

.261 

.115 

—.008 

—.030 

—.016 

.004 

.002 

.001 

—.001 


.284 

.239 

.136 

.035 

—.019 

—.026 

—.011 

.000 

.003 

.001 


r 
P 


.707 
.229 


.645 
.122 


.588 
.0622 


.558 
.0411 


.508 
.0195 


.707 
.446 


.625 
.239 


.563 
.0829 


.532 
.0365 


.486 
.0134 



It will be seen that although the coefficients I of the two original formu- 
las (11) and (12) differ from each other very considerably, yet after 12 ap- 
plications they approach pretty nearly to each other and also to those of the 
12-fold formula of Table I. The idea is thus suggested, that when any 5- 
term formulas similar to (8), (11) and (12) are applied a large number of 
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times in succession, the resultant coefficients approximate to a common form, 
and that they resemble the coefficients of Table B. a good deal more closely 
than those of Tables A. or C. This may perhaps be regarded as an argu- 
ment in favor of the use of the formulas of Table B. for purposes of adjust- 
ment. It is reinforced by the fact that when other formulas are repeatedly 
used, they take a similar outline. The 9-term formula of Table A., in ex- 
act numbers, is 

M 0=2-Tr[ 59w 0+ 54 ( M l+ M -l) + 39 ( W 2+ M -2) + 14 ( W 3+M- 3 )— 21(M 4 +M_ 4 )], 

and if it is used four times in succession, the coefficients of the equivalent 
formulas are, to three places of decimals, 

.202 .185 .141 .086 .035 —.001 —.018 
—.018 —.011 —.004 .001 .002 .001. 
These differ but little from the coefficients in the 17-term formula of Table 
B. Again, the 13-term formula of Table A., in exact numbers, is 

< = tM 25m o+ 24 K+ m -i)+21K+w_ 2 )+16(w 3 +«_ 3 ) 
+ 9(m 4 +w_ 4 )+ 0(u 6 +«_ 6 )-ll(u 6 +«_,)], 

and, after four successive applications of it, the resultant formula is found 
to have the following coefficients : 

.138 .132 .117 .095 .069 .043 .021 
.004 —.008 —.013 —.013 —.010 —.006 
—.003 —.001 .001 .001 .001 .001. 
There is no very great difference between these and the coefficients of the 
25-term formula of Table B. 

A comparison of the values of r and p in Tables I. and II. affords some 
rather unexpected results. Although the original formula (12) gives values 
to r and p as large as, or larger than, either of the original formulas (8) and 
(11) do, yet after twelve successive applications, it gives smaller values than 
the others do. In other words, formula (8), which was found by the method 
of least squares, and makes, when the conditions assumed are fulfilled, a 
more accurate single adjustment than any other 5-term formula, is not the 
one which gives the greatest accuracy when the adjustment is repeated sev- 
eral times, or even when it is repeated once only. Likewise formula (11), 
which under the same conditions makes the smoothest single adjustment of 
any, is not the one which gives greatest smoothness when the adjustment is 
repeated. 

Other calculations have been made for the purpose of learning what the 
results will be when the 25-term formulas of Tables A., B. and C. are used 
in succession, that is, when any one of them is used twice, or when any two 
are used one after the other. Each of the complete resultant formulas in- 
cludes 49 terms, with 25 different coefficients, and from 17 to 25 of these 
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are large enough to affect the third decimal figure. It is thought unnecessary 
to give any of them here in full, but the ratios r and p for all of them are 
shown in Table III. For instance, when the 25-term formulas of Tables 
A. and B. are used in succession, no matter which is used first, the error- 
ratio and irregularity ratio of the resulting adjustment will be 
r = .281 p = .000191. 

Table III. 





"Values ol 


r. 


Values of 


p. 




A 


B 


C 


A 


B 


C 


A 


.266 






.002649 






B 


.281 


.313 




.000191 


.000324 




C 


.272 


.291 


.279 


.000301 


.000185 


.000313 



When the 25-term formula of Table C. is applied twice, the ratios for the 
result will be r = .279 p = .000313, 

and so on. As all the values of r in Table III. are larger than .250, it will 
obviously be difficult, if not impossible, ever to adjust a series by these 
methods so as to reduce the error-ratio below J. It diminishes but very 
slowly with each repetition, the ratio for the second application of A., for 
instance, being .266, while for the first application Table A. shows it to be 
.300, and we have 

.266 __ eQ7 

so that the probable error of the second adjustment is little less than nine 
tenths that of the first one. 

We also see that the smallest value of p in Table III. is not, as might 
have been expected, that which results from two applications of the B for- 
mula. Smaller values are given not only by two applications of the C, but 
also by either A and B, A and C, or B and C. 

Signs of Eesiduals.— In the Analyst for Jan. 1878, p. 3, it was 
shown by me that when a series has been adjusted by a formula whose error- 
ratio r is small, a safe limit being r < \, the probability that the signs of 
any two adjacent residual errors will be alike is approximately 

1— 2r 2 



and as a consequence of this, that the most probable number M of signs fall- 
ing in groups of only one or two signs each, in a periodic series, is given by 

Jf=(p 1 +2p 2 )^±.6745 1 /[p 1 (l_p 1 )+4p 2 (l- i32 )] v /JV, } ( 14 ) 
where iVis the whole number of terms in the series. The probable num- 
ber falling in groups of more than two is 

W = N—M. 
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The above will answer also for non-periodic series, if JV-§ is put in place 
of N, in which case the first and last signs of the series, and as many adja- 
cent ones as happen to be like them, are not counted. If we write 

M' = aN±bvN, 
the values of a and b, for different values of r, will be as here given in 
Table IV. 

Table IV. 



r 


a 


5 


a 1 


h 


.25 


.450 


.55 


.467 


.34 


.30 


.426 


.56 


.451 


.34 


.35 


.396 


.56 


.430 


.33 


.40 


.359 


.57 


.405 


.33 


.45 


.314 


.58 


.373 


.33 


.50 


.259 


.58 


.333 


.32 



The process by which (14) was obtained assumes that any four consecu- 
tive signs may be regarded as independent of each other, except as adjacent 
ones are connected by the relation (13). This is only approximately true. 
For the purpose of testing an adjustment, it may be better, and is simpler, 
to consider only the number of permanences of sign, that is, the number of 
cases where adjacent signs are alike, in passing along the series of residuals 
from one end to the other. The most probable number of permanences will 
be, in a periodic series, 

P = qN±£745i/[q(l—q)Nl (15) 

and the same for a non-periodic series, only putting N — 1 in the place of N. 

If (15) is written 

P = a 1 N±b 1 - l /N 
the values of a : and b 1} for different values of r, will be found in Table IV. 
For example, in the observed series of signs at p. 5 of the Analyst refer- 
red to, we have r = .372, N = 70 and by interpolation in Table IV., 
a = .380, b = .57, a x = .419, b x = .33, so that computation gives 
M' = .380 x 67.5±.56 1 /(67.5) = 25.7±4.6, 
P =.419x69 ±.33]/(69 ) — 28.9±2.7. 
Observation by actual count in the series shows 

M' = 22, P = 25. 

Both fall short of the computed numbers, but the deviation is on the safe 
side, and must be considered accidental. It is only when the observed val- 
ues of M! and P exceed the computed values, by more than the probable errors, 
that we are to infer the existence of a danger that the adjusted series has been 
smoothed out too much, so as to deviate materially from its true law. 



